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SUMMARY 


A solution is obtained to the transient- temperature distribution 
in a semi-infinite two-component composite slab of arbitrary materials 
subjected to an instantaneous application of aerodynamic heating with 
constant equilibrium temperature and heat-transfer coefficient. The 
numerical results are tabulated in a form to permit easy computation of 
heat-transfer problems typical of aerodynamic testing. The solutions 
are valid for finite two-component slabs as long as the times considered 
are small compared with the diffusion time of the backing material. 

Analytical results obtained from these solutions can be used to 
determine (a) the heat-transfer testing time for ’which the outer skin 
may be assumed to act as a calorimeter without exceeding a given error 
or (b) correction curves by which the Indicated calorimeter heat-transfer 
coefficient may be multiplied to obtain the true heat-transfer coeffi- 
cient. For such a correction curve to be valid, the bond between the 
two materials must have negligible thermal resistance, a condition dif- 
ficult to attain if the slab is not composed of two metals. 

Since the differential equation for the temperature distribution is 
linear, the principle of superposition is valid. Consequently, the prob- 
lem of continuously varying equilibrium temperatures and heat-transfer 
coefficients may be treated by using the tabulated solutions and con- 
sidering the continuous variation as a series of superimposed step 
functions . 
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INTRODUCTION 


The increase in flight speeds for "both operational and experimental, 
missiles and airplanes has resulted in renewed interest in the problems 
involved in fluid- and solid-state heat conduction. The two problems of 
transient behavior of a two -component slab, consisting of an outer skin 
with a backing material, exposed to aerodynamic heating that have become 
of prime importance are: (l) the performance of the surface as measured 

by the ability of the material to withstand extreme heating rates on the 
full-scale vehicle during the course of its mission and (2) the transient 
testing techniques presently in wide use wherein the aerodynamic heating 
is varied and the resultant time -wise variation of the temperature of a 
skin, either with or without a backing, is then used to obtain heat- 
transfer rates and coefficients. For example, these transient testing 
techniques are used for shock tubes, free-f light pilotless rocket- 
propelled models, the sudden insertion of a model Into a wind tunnel, or 
the rapid variation of the stagnation conditions in wind tunnels. 

In the interest of simplicity, these techniques, except that used 
for the shock tube, customarily assume that the outer skin of the slab 
acts as a calorimeter and absorbs all the aerodynamic heat with negli- 
gible loss to the backing material. However, even when the backing 
material has very low conductivity, heat is continually being transferred 
across the interface and this heat transfer may introduce large errors in 
the answers obtained by the calorimeter assumption. 

Two recent papers (refs. 1 and 2) treat the heat transfer to such a 
semi-infinite composite slab. Reference 1 is applicable to shock- tube 
testing wherein a very thin metal plating is formed on an insulator and 
the temperature response of the metal is used to determine the heat 
transfer. This study Is applicable to arbitrary heat flux rates -but 
assumes the plating to have such negligible thickness that the interface- 
temperature solution may be found as a perturbation to the surface solu- 
tion in the absence of the plating. The solutions obtained are valid for 

times much greater than the diffusion time in the plating."*" Reference 2 
considers the case of the heating of a composite slab in which the heat 
transfer is proportional to the difference between wall and equilibrium 


In the discussion of transient heating problems It is convenient to 
have reference values of time such as "diffusion time" and "relaxation 
time" for comparison purposes. The diffusion time is the quotient of the 
square of a typical length (for example, the thickness of a slab) divided 
by the thermal diffusivity of the material. The relaxation time, appli- 
cable in exponential decay behavior, is also called the time constant and 
is that value of time which makes the absolute value of the exponent of 
e unity. 
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temperatures. The equilibrium temperature is a step function. The outer 
skin has finite thickness but infinite conductivity; thus, the temperature 
in the skin is uniform at any given time. A solution is presented for the 
heat conduction to the backing material, no coupling between this loss of 
heat to the backing and the temperature distribution in the outer skin 
being assumed. This solution is valid for times that are small compared 
with the relaxation time of the outer skin. Reference 2 shows that con- 
duction even to relatively good insulators introduces appreciable error. 

The present investigation was initiated to obtain solutions for the 
heat transfer to the composite slab without these restrictions. The mag- 
nitude of the errors arising from the calorimetric assumption and pos- 
sible correction factors were to be determined. When the thickness of the 
backing material is such that the diffusion time for the backing material 
is much larger than any of the testing times, the backing material then 
acts as if it were of infinite extent during the test. Consequently, solu- 
tions to the problem of a semi-infinite composite slab would apply under 
this condition. The problem considered in this report is the transient 
temperature distribution in a semi-infinite two-component slab subjected 
to aerodynamic heating at a rate equal to the product of a constant heat- 
transfer coefficient and the difference between surface and equilibrium 
temperatures . The equilibrium temperature is assumed to be a step function 
and the outer skin is assumed to have both finite conductivity and thick- 
ness. Solutions are found which are applicable for times varying from 
much less than to much greater than the skin diffusion time. The super- 
position principle permits extension of the results to arbitrary timewise 
variation in heat-transfer coefficient and equilibrium temperature. The 
investigation reported herein was conducted in the Gas Dynamics branch of 
the Langley Aeronautical Laboratory. 


SYMBOLS 


a,b 

A,B,C,D 

c 

H = h/k 
h 

A 

h 


constants 
constants 
specific heat 

functions defined by equation (23) 

heat-transfer coefficient 
indicated heat -transfer coefficient 

(dT/dt) n 

indicated heat-transfer coefficient at | = 0; pci - 

3 7 m m 
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A 

h l 

indicated heat- transfer coefficient at 

H 

11 

UJl 


- 

J 

coefficient of g in expression for I 

? and equal to +1 


* 


or -1 




k 

thermal conductivity 




T lh 
K 1 

lH l 




Ln 

Laguerre polynomials (see appendix A) 




l 

distance from surface to interface 




m,n 

integers 




P 

parameter in inverse LaPlace transform 




P 

variable of the Laplace transform 




Q 

heat flux to surface per unit area per 

unit time 



q = ^i/a 




• - 

s = p/ka 2 





T 

temperature 



• 

*e 

equilibrium temperature 




T o 

temperature at g = 0 




T 

relax 

relaxation temperature 


9 


"^relax 

relaxation time 




t,t' ,u 

time variables 



.... 

X 

coordinate measured normal to surface exposed to heating 



a 

thermal diffusivity, k/pc 




p = 1 - 0 





1+0 




- 
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6 = 0 - 2.5 

F solutions to P tan r = L 

s s s 

x = E to 

t = x/Z 

p density of material 

a = = J P 2 C 2*2 

tj , X, r dummy var iable s 

} = t + Leu 

t - 2. 
y 2co 

O p 

n - jr - + 

cut 
CD = I q ' 

Z 2 


inverse Laplace transform 

A bar over a symbol denotes the Laplace transform. Except as desig- 
nated in the preceding list, the subscripts 1 or 2 refer to evaluation in 
region 1 or 2. 



THEORY 


General Solution 

The problem to be solved is the one -dimensional transient tempera- 
ture distribution in a semi-infinite two-component slab which has one 
face suddenly exposed to a fluid. This fluid will transfer heat to the 
exposed face at a rate proportional to the difference between the constant 
equilibrium temperature of the fluid and the outer surface temperature of 
the slab. In the interval 0 < x < Z, the slab is composed of a material 
with properties p-^, c^, k l> and and, in the interval Z < x < «>, 

it is composed of a material with properties Pg, Cg, kg, and ag. No 
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thermal resistance is assumed at the interface in the slab. The mathe- 
matical formulation of this problem is given in terms of the following 
governing differential equations and boundary conditions (see fig. 1): 


o T-, t 5Ti . . . 

"° (t>0 ' 0<x< (1) 


8% 

8x 2 



(t > 0, l < x < <») (2) 


^ = t 2 = 0 


Tg ** 0 


8t 


k l 8T + h ( T e “ T l) - 0 


rp _ qi 

J T 2 


. dTp 

k isr’ k 2^ s 


Sx 


(t < 0, 

all x) 

(5) 

(t > 0, 

x <») 

00 

0 

It 

t > 0) 

( 5 ) 

(x = l, 

t > 0) 

(6) 

(x = l, 

t > 0) 

(7) 


The LaPlace transforms T (p) = / e”^T (t)dt and the definitions 

Jo 

of q s Vp/a are then employed to reduce equations (l) to ( 7 ) to the fol- 
lowing forms : 


d 2 ^ 


- 1 ± \ = 0 


(0 < x < l) (la) 


d 2 T, 


— 2^ - = 0 

dxr 2 2 


(l < x < w) (2a) 
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-»0 

(x -» eo) 

(te) 

, dTn - ML 

k n — i - hT n + — £ = o 
- 1 - dx 1 P 

(x = 0) 

(5a) 

H*' 

II 

ro 31 

(x = l) 

(6a) 

dT n dTo 

k. — i = ko — £ 
1 dx ^ dx 

(x = 1 ) 

(7a) 

The general solutions to equations 1(a) and 2(a) 

are: 


T x = Ae-^ + Be^* 


(8) 

\ = Ce^ + De^ 


(9) 

The value of D must he identically zero to satisfy equation A (a) , 
Application of the boundary conditions (eqs. 6(a) and 7(a)) yields 

♦ 

B = pAe“ 2q -l Z 


(10) 


C = (1 + p)Ae~ 9j - l+q21 (11) 

The remaining constant A is then found from equation 5(a) to he 


H,T -j 

A = 1 e 1 _ 

P q.! + H] 


1 - , e ' 2 *! 1 

V ii + % 


-1 


( 12 ) 


The solutions to the differential equations in the transformed plane 
become 



Ei i 1 [e"^l x + 

P ^1 + H 1 x ll - Hi ^-2q-| l L J 

*1 + Hi 


(0 < x < l) 


( 13 ) 
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Hr 

F 


% 1 

P q-L + Hj 


1 - B 


1 H- P 

*1 - H x -2a 


“4l 


L.* i-p fx 7\" 
* + k£ I+p (x " l) _ 


(x > l) 

m 


y’-i/T-A 

The inverse LaPlace transform/ \7jf^J of equation (13) is known 


for the special case a = 0, P = 1 which corresponds to a perfect insu- 
lator as the backing material. If the exponentials of equation (13) are 
expressed in hyperbolic form, the inversion on page 259 of reference 3 
applies . 



= 2L 


00 

I 


2 2 

. + IT + l 

s=l S 


cos r s (i - 6 )1 
cos r 




,1 - e 


IF 


(15) 


where P„ are the roots of the equation 


T tan r = L (16) 

S 8 

If the slab is composed entirely of the same material ^p = 0$ 
k^ = kg), equations (13) and (l4) become identical and have as their 
transform 


-S- = erfc 



Hqx+aHq^t 
e erfc 



+ 



(17) 


For the general case of a ^ 0, P ^ 1, the inverse transform of 
T T as expressed in the form of equations (13) and (l^) is neither 

known nor easily determined. Recourse is then made to expansion of the 
right-hand side of these equations. (Rote that |p| ^ 1.) An expansion 
of this type generally results in solutions which will converge more 
rapidly for small values of time. 
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Since 


-1 


e " 2 ^ ) „ V [p e ~^ 

q i + H l 1 ^\*1 + E 1 


- H-, “2q. n l 


,1 - 0 


(18) 


substitution of equation (l8) into equations (13) and (ill-) yields 

Ti p n (q.i - H x ) n _-'l 1 U2n+|) p n+l _ Hl ) n _-a 1 l[ 2 (n+l)-(] 

tT “ H ! A ~ 7 „ N n+l C + P 7 7 \n+i e 

6 n=0 ^1 + H l) 'll + H l) 


h u (1 , B) f ^ (*1 + 

^ - V 1 + s) 1 T TTTTpi e L J 

n=0 v 1 x/ 


(19) 


(20) 


If the orders of summation and integration of the inversion process 
are assumed to be Interchangeable, then 


- 1 /- 


00 _1 




n=0 


f (ii - %) n - v (a, * 6) . 

p 7 TEb I e + 

(ii + E i) 


g a+1 fa - S ± f ^l^Cn+D-s] 
P + 


, n+1 


( 21 ) 


r 


-1 't 


CO -1 




n=0 


; fa - gif e -v[ 2n+1+ | w"- 1 ’] . 

P , _ sn+l 


fa + H i) 


( 22 ) 


Let the function G n (P,cu,L) 


be defined as follows: 
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o O 00 v 2 p 2oIj(x - ) _rj 

G (P,o>,L) = -p / e" / V 2 “4 \(2i})dr) dX (23) 
\pw p/2co JO 

In appendix A it is shown that the inversion of the general term of 
the series has the form 


= G n (P,o>,L) 


% (q.1 - Hx) 11 

P / TT \*1+1 

(»1 + H l) 

Substitution of equation (24) into equations (21) and (22) produces 


X 


(24) 


^ = Y p n G n (2n+l ^ L ) + 0 n+ ’'"G n { 2 ( n+1 ) - £ lJ 


n=0 


(25) 


09 

^ - (1 + 0) £ p\ 2n+l + ^ i=|(g-l),cu,L 
n=0 


(26) 


The formal integrations required for the evaluation of G n for a 
range of values of n from 0 to 5 are given in appendix B. The 

/ dG n aGn\ 

associated timewise and spacewise derivatives — ry and are also 

\8co dg / 

presented. These derivatives permit the evaluation of both the correc- 
tion curves to the indicated calorimetric heat-transfer coefficient and 
the amount of heat crossing a given station, in order that estimate of 
the flux through the backing material may be obtained-. 


Significance of parameters L, X, and 3 
The physical significance of the nondimens ional parameters L, X, 

and p will be briefly discussed at this point. The parameter L = — 

ki 

is proportional to the ratio of the temperature difference required to 
transport a given amount "of heat per unit area across the slab to the 
temperature difference required to transport the same amount of heat per 
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unit area from the fluid to the slab . This may "be illustrated for the 
steady-flow case as follows : 

Q = h[l e - T(0)J = T iQ) . - ■ T ( x ) . 


Therefore , 


T(0) - T(l) = hi = L 

T e - T(0) ^ 

Now if L is small, a small temperature difference across the outer 
skin is indicated; thus, the surface temperature and heat transfer when 
T(o) -»T e will he sensitive to changes in the temperature at the inter- 
face. Since the interface temperature is in turn influenced strongly by 
the backing material, the relative influence of the backing material 
increases with a decrease in L ' for times of the order of the relaxation 
time. (For very large times the backing material must always be the 
dominating factor.) The significance of L may be realized in an alter- 
nate way by considering that, when L is small owing to a thin skin, high 
thermal conductivity, or a low heat-transfer coefficient, a greater por- 
tion of the heat transferred to the surface is available for transfer to 
the backing material and the backing material can exert a larger influ- 
ence on the temperature time history of the skin for finite times. A 
plot of the parameter L obtained for various skin thickness and material 
combinations with h = 0.1 Btu/sq ft-sec-°F is shown in figure 2. A range 
of L from 0.001 to 0.5 appears to cover practical outer skins which 
might be used for high-speed vehicles . 

The significance of the parameter X = £ cd =,p£— is most easily 

V pel 

described from consideration of the aerodynamic heating of a slab of 
thickness l having infinite conductivity (8 t/8x = 0) . The slab acts 

as a calorimeter absorbing all the heat input. This slab is initially 
at T = 0 and is subjected to a step input in equilibrium temperature at 
t = 0. The temperature of the slab at time t may be expressed as 



Thus, the value X = 1 corresponds to the relaxation time in such a 
slab calorimeter. 

For a perfectly insulated skin having finite conductivity and thick- 
ness, the temperature distribution departs from that of the calorimeter 
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as L increases from zero. However, at the value A = 1, the tempera- 
ture is approximately the relaxation value given by T relax = T e (l - ij. 

Since t « A , the ratio of any time t to the relaxation time t relax 
is simply A^. 

The parameter (3 denotes the relative insulating perfection of the 
backing material. A value of (3 = 1 applies to a perfect insulator j 
P = 0 describes the homogenous composite slab constructed of the same 
material throughout; and |3 = -1 describes a perfect conductor which 
maintains a constant temperature at the interface at | = 1. Values of 
0 appropriate to various combinations of outer skin and backing material 
are shown in figure 3- The outer skin materials have common abscissas 
and the intersection of these vertical lines with the curves for the 
backing materials determine the appropriate value of 0. If the materials 
for skin and backing are interchanged, the sign of 0 is reversed. 


Particular Solutions 

Since many aerodynamic heat-transfer experiments are performed with 
a model which has a skin of high conductivity and density compared with 
that of the backing material, it is desirable to have a form of solution 
applicable to these circumstances . The solution for a perfect insulator 
backing material is either the series given in equation ( 15 ) or that 
given in equation (25) when a value of 0=1 is used. These alternate 
solutions converge most rapidly at opposite times . The number of terms 
required to express the timewise derivative for equation ( 15 ) is large 
as t -» 0 and decreases rapidly with an increase in time; thus for values 
of t greater than about one-half the relaxation time, only one or two 
terms are required to give the value of T-j_ jT e and its timewise and 

spacewise derivatives with acceptable accuracy. Conversely, although 
one term is sufficient for equation ( 25 ) when t -» 0 , the number of terms 
required increases with time. A series thatnilll converge rapidly for 
the times of interest for aerodynamic testing purpose’s when 0 1 may 

be obtained by applying the solution of equation ( 25 ) as a perturbation 
to the solution of equation (15). If T-^Tg for 0 = 1.0 is added and 

subtracted to the right-hand side of equation ( 25 ), the resulting equa- 
tion is 

00 ( -x 

% = (l ) f . 1 ' Z r" » n V aM -S'‘"' L > + V- ■ 0 n+1 K[2(n+l)-5,«>,L] ■ 

n =0 ' 


( 27 ) 
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where is evaluated from equation ( 15 ) for the appropriate 

' e / p=l 

values of L, ca, and g . 


Of particular interest is the temperature history on the surface 
exposed to the transferring fluid (g = 0) and that on the interface (g = 1) 
between the two different materials. Substitution of g = 0 into equa- 
tion ( 27 ) produces the following set of equations: 


*L( 


5=0,a),L) 



£ (l - p n )G n (2n ) m J L) + £ (l - 
[n=0 n=l 


( 28 ) 



A corresponding set of equations may be obtained for the interface 
by substituting g = 1 into equation ( 27 ) : 



The evaluation of the terms of equations (25) and ( 26 ) becomes 
progressively more difficult as n Increases. Consequently, if the 
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temperature at the interface (| = 1) has to be evaluated for any reason 
or if the values tabulated in this report are to be used, then the tem- 
perature distribution in the backing material can advantageously be 
obtained by Duhamel's method. 


Since the temperature distribution in a semi-infinite slab (l < x < °o) 
when one surface at x = l undergoes a unit- step-function increase in tem- 
perature is (ref. 3) 

T (x,t) = erfc (34) 

2f^t 

Duhamel's technique yields 

(x-Z) 2 

T 2( (x-Z)(l+jl) P t e ’4 <X 2 (t-r) 

* 2f~ Jo (t - t) 3/ 2 


X ?n °n( 2n+1 ’¥' L ) 

n=0 \ i / 


dT 


(35) 


RESULTS AND DISCUSSION 


General Results 


The values of G n and its derivatives for n = 0 to 5 and 
|=0 and 1 were computed on a card-programed calculator. Details of 
this procedure are discussed in appendix C. Hie results of these com- 
putations are given in table I for values of L between 0.001 and 0.5 
with 0.02 A = 8. Table I is arranged for future computing convenience. 
Values on the left for | = 0 are alined so that pairs in G or its 
derivatives on the same horizontal line are operated on by the same power 
of p as prescribed by equations (28) to (30). 


As discussed In appendix C the computing procedure at times involved 
the subtraction of two large and nearly equal numbers obtained by approxi- 
mations. Certain voids appear in table I where an inspection of the 
machine results indicated insufficient accuracy. There was no rigorous 
rule for determining which answers to discard, but all answers with less 
than four significant figures were usually eliminated. Those with less 
than three significant figures were always eliminated. 

Typical values of G n and its associated derivatives which are 

required at | = 0 are plotted, against L for A = 1.0. (See fig. 4.) 
The values shown are those needed in. the evaluation of the first series 
in the braces of equations (28) to (30). It is immediately obvious that 
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the convergence of the series is much more rapid for large values of L 
than for small values of L. For example, at L = 0.5 the ratio of 
to Gq is of the order of 10“6 whereas at L = 0.001 the ratio is of 
the order of one. 

An indication of the accuracy which results from consideration of a 
finite number of terms to represent the infinite series of equation ( 25 ) 
may be obtained from figure 5- In this figure are plotted the values of 
the ratio of T]_/T e for (3 = 1 when the upper limit of s umma tion in 

equation ( 25 ) is reduced from infinity to 0, 1, 2, 3* 4, or 5 to that 

obtained from table H . ( Table II contains the values of Tq/T e as 

determined from equation ( 15 ) by using a sufficient number of terms to 
get the desired accuracy.) The ratio for values of L = 0.001, 0.1, 
and 0.5 is plotted as a function of A. Hie number of terms required 
to reduce the error to a given limit is seen to increase with an increase 

in A at a constant L and to increase with a decrease in L at a con- 

stant A. For a constant L, the behavior of T-j/T e for any finite num- 
ber of terms becomes oscillatory as A increases. The errors for 
0 p < 1 at a given A and L should be less than those indicated in 

figure 5 for P = 1 since the terms neglected contain p n or p 
In other words, the neglected terms for 0 P < 1 have a smaller con- 
tribution than those for p = 1.0. 


Illustrative Example 

The effects of the backing material in a composite slab under typical 
transient test conditions will be illustrated by the following study. ■ The 
numerical results obtained are found by using the values of table I for 
0 Is n % 5 and those of table II. Hie mathematical assumption of an infi- 
nite extent of backing material to represent the finite backing of the test 
is valid for values of p near unity as long as the times considered are 
much less than the diffusion time of the backing. The diffusion times, 
which are equal to the square of the thickness of the backing material 
divided by the diffusivity of the backing material, for l/k- inch.- thick 
balsa and mahogany are approximately 120 and 170 seconds, respectively. 
Consider two stainless-steel outer skins (k = 0.0028 Btu/sec-ft-°F, 
a = 5*2 x 10-5 sq ft/sec) with thicknesses of 0.060 inch and 0.030 inch 
backed by either a perfect insulator (p = 1 . 000 ), balsa (p 0 . 975 )* or 
mahogany (p « 0.950)* These composite slabs are subjected to a step 
input in equilibrium temperature (for example, by sudden exposure to the 
airstream of a wind tunnel) and the heat-transfer coefficient h is 
O.O 56 Btu/sq ft-sec-°F. The values of L are then 0.10 and 0.05 for the 
0 . 060 - and 0 . 030 -inch thicknesses, respectively. 
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The resulting behavior of the outer skin is illustrated in figure 6 
which shows the temperature distribution at the outer face ( 4 =0) as a 
function of time. The effect of the heat lost to the backing material is 
evident in a lower temperature in the outer skin at a given time, a larger 
effect being observed for the thinner skin. 

An indicated heat-transfer coefficient is often obtained from the 
temperature-time history of such an experiment. This heat-transfer coef- 
ficient is determined by assuming that the slab acts as a calorimeter with 
an infinite value of k; that is, 

T e - i 

A 

Since the value of h depends on the value of 4 at which T is meas- 
ured for finite k, the term h Q or h^ will be used to designate the 

values of h which would be found when T and dT/dt are evaluated at 
| = 0 or 1. 

The ratio of the indicated heat-transfer coefficient to the true 
heat-transfer coefficient h/h is plotted as a function of time in fig- 
ures 7 and 8* Figure 7 shows the ratio at 4=0 and 1.0 for the 
0.030- inch- thick skin. Except for the initial discrepancy at early times 

^ A A 

the values of h Q and h^ are nearly equal. At t = 0, hQ is infinite 
while h-, is zero, but after about 0.1 second the difference is minor. 
Figure 8 shows the ratio hQ/h as a function of time for both the 
0.030-inch and 0.060- inch skin. It is evident that h Q = h only at a 

single discrete point for each combination of skin thickness and backing 
ma terial. Furthermore, even the perfect insulator does not give a ratio 

of unity for h/h except at one value of time. It may be shown for the 
perfect insulator that the asymptotic value of h/h is 


!im £ 
t -» » 



L 


1 - 


1 T 

5 L 


1 + |L 
P 


1 + 


II 


1 1 2 


» 1 


1 

3 


L 


(L « 1) 
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Thus even the perfect insulator gives errors of 5/3 and 10/3 percent 
for the cases considered (0. 030-inch and 0.060-inch-thick skins). The 
large increase with t ime of the deviation of h/h from unity for the 
imp erfect insulators is striking. This error increases with time because, 
as the heat transferred into the skin at 1=0 decreases monotonically 
with t ime , the heat transferred out at | = 1 first rises with time to a 
rrw. vi mi im (at approximately A = 1.0) and then subsequently decreases but 
at a rate slower than that at | = 0. An approximate method for esti- 
mating the amo unt of heat conducted across the interface for the condi- 
tions of A « 1 and k -» co was also presented in reference 2. 

If an error of 10 percent were set on the deviation of indicated to 
true heat-transfer coefficient, the time limit for useful data with the 
respective backings composed of mahogany and balsa would be about 1 sec- 
ond and 2.5 seconds when 1 = 0.030 inch and 2 .^ seconds and 5*7 seconds 
when l = 0.060 inch. If, instead of using a time limit, a temperature 
rise limit for a 10 -percent error were desired, it may be obtained from 
a plot similar to that of figure 9 . For l = 0.030 inch, temperature 
rises greater than 30 percent of T e for a mahogany backing or 65 percent 
of T e for a balsa backing would give values of h/h deviating from 

unity by more than 10 percent. Approximately the same limits also apply 
for l = 0.060 inch. 


CONCLUDING REMARKS 


A solution has been obtained to the transient temperature distribu- 
tion in a semi-infinite two-component composite slab of arbitrary mate- 
rials subjected to an instantaneous application of aerodynamic heating . 
with constant equilibrium temperature and heat- transfer coefficient. The 
numerical results are tabulated in a form to permit easy computation of 
heat- transfer problems typical of aerodynamic testing. The solutions are 
valid for finite two-component slabs as long as the times considered are 
smal 1 compared with the diffusion time of the backing material . 

Analytical results obtained from these solutions can be used to 
determine (a) the heat-transfer testing time for which the outer skin may 
be assumed to act as a calorimeter without exceeding a given error or 
(b) correction curves by which the indicated calorimeter heat-transfer 
coefficient may be multiplied to obtain the true heat-transfer coefficient. 
For such a correction curve to be valid, the bond between the two mate- 
rials mus t have negligible thermal resistance, a condition difficult to 
attain if the slab is not composed of two metals. 

Since the differential equation for the temperature distribution is 
linear, the principle of superposition is valid. Consequently, the problem 
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of continuously varying equilibrium temperatures and heat-transfer coef- 
ficients may he treated by using the tabulated solutions and considering 
the continuous variation as a series of superimposed step functions . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June k, 1958. 
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APPENDIX A 

INVERSION OF SEE LAPIACE TRANSFORM 


The steps required to invert the general terms of the series of 
equations (21) and (22) will be determined in this section. The inver- 



(Al) 


will first be found and then the constants a and b will be assigned 
pertinent values. 

Let s = -2—. It is then necessary to find the inversion for 
ba d 


-1 

X 



-1 

- y 9 

-2ab>/s 

1 1 e 

(afr - D n " 

X 2b + D n+1 _ 


4a 2 ]fs 

(2\Jb + l)”* 1 


(A2) 


The square-root image relation (see page 123 of ref. 4) can be used to 
evaluate equation (A2) . Thus if 

- 1 . 


f(t,) 'X t g(s) ] 


then 


co u -1 

:== f e lrt ' f(u)du = -p g(Vs) 

V«t ' JO \/s 


(A3) 


and the evaluation of the following expression is required: 

- 1 . 




e -2abs (2s - l) 

8 (2s -K l) n+1 


n “ 

+1 
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03ie term 


and 



(0 < t < gab) 

► 

(2ab < t) 

- 


(A4) 


-ir 


X 


(2s - l) n 


(2a + 1) 


n+1 



(A5) 


(See page 298 of ref. 5 and page 129 of ref. 4.) 03ie symbol L^t) 
designates Laguerre polynomials of the following form: 


L n (T) 



■h 




L n ( -> = I At If 

m=0 


m 


(A 6) 


Application of the convolution theorem to equations (a 4) and (A5) results 
In: 



/ S n 

(2s - 1) 
(2s + 1) 


' U " 2ab -t/2 t 


\ e ~ T/ tl 4 n ( T ) (l)^ T + 

p U _ t / 2 

5 / e L^t) (0)dr 

d J u-2ab 


P u-2ab _ T /p 

1 I e ‘ L (r)dr 

2 Jo 


(AT) 
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Combination of equation (A7) and. the image formula (A3) yields 

u2 




_i_ 

fs 


-2abv/s 


(2x/s - 1) J 


(2v/? + l) n+1 


= _1_ f ' f 

\IH t 7 J 2ab J 0 


u-2ab -t/2 

L n (f)dT du 


(A8) 


For s = cp, ct’ = t, it may be shown that, if 

-1 


„.,r [s( s )]= f ( t! ) 


then 


-1 


at [ s(pc) ] - I f (£) 


(A9) 


Consequently, since c = — i-, equation (A8) may be operated on to yield: 

4a 2 


2T 


" 

1 

—L 

e 


pE 

. ika. 2 ) 

a " 


\ 

4a c 

( 2 

J5* 1 ) 

- 


ka 2 


l jfs4a^t 


335 r-** e -\ (T)dT du 


r e 's^t r 

J 2ab u 0 


or 


-1 




- b /p 


(v/p - a) 


n 


(/p + a) 


n+1 


If . a = H Ja and b = equation (A10) becomes 
“ '/a. 


— f 

4a /ict J 2ab 


u 11 - 

00 "l6a^t o u-2ab -J 

e / e ^(TjdT du 

0 


(A10) 
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-IPq 


( q 


( q + H) 


H ) 

aFI 


1 — r 

i( 'H 1 Jjtat J ZISF 


dt r- 2UIP e- T /%( T )d T du 

J 0 


(All) 


If 


and 


n = t/2 



u 

EiSj 


are substituted in the right-hand side of equation (All ) } the following 
equations result: 



e - lF ^ fa- H i) n 

p / _ \n+l 

+ H l) 


_ 2 _ 





L n (2t 1 )dT 1 dX 


= G n (P,a),L) 


(A12) 
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APPENDIX B 


EXPRESSIONS FOR THE VARIOUS G_, — AND — 2. 

n &> 2 d| 

FOR VAIUES OF n FROM 0 TO 5 


Formal integration of the equations for the G n for values of 

from 0 to 5 and subsequent partial differentiation with respect to co 
and g produced the following expressions: 


G 0 = erfc \|r - e^erfc 0 
= -JLe^erfc 0 

dG 0 L 3 / 2 -t 2 t 2 ft _ a 

acu 2 ^ A 

G = -erfc i*r + (l - 4 A\/l 0) e^erf c 0 + — A\/Le 
x \fn 


-* 2 


Sg 

aT 

Sg. 


-V 2 


= -Lj fl + AA\/L0^e^erfc 0 + j -p AL^/ 2 e 
7 V* 

^ = -L 2 ( 5 + WE0> fl e rfc 0 + + ^e^ 

G 2 = erfc t - (l + 8A 2 I0 2 + 4A 2 Lje 2 erfc 0 + A 2 L0je _ ^ 2 

= -Lj(l + 8A\/I0 + Q-k 2 !# 2 + 4A 2 lJ e 2 erf c 0 + M AL 3 / 2 (l + A^0)e“^ 
= -L 2 ^5 + 16A0L0 + 8A 2 L0 2 + ^A^Je^erfc 0 + 


_L_ 

& 


^ + 8AL 3 / 2 (2 + A>/L0) 


-* 2 


ro p 
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= -erfc \|f - e 2 erfc 


0 ^-1 + 4a/l0 + 8A 2 L0 2 + 4 A 2 !, + y. A 3 L 3 ^ 2 0 3 + 
l6A 5 L 5//2 0) + ^=(Wl + 8A 2 L0 + y A 5 L 5 / 2 0 2 + ^ A 3 ! 5 / 2 ^"* 2 


^ = -e Q erfc 0 j (l + 12AL 3 / 2 0 + 24A 2 L 2 0 2 + 12A 2 L 2 + ^ A 3 L 5 / 2 0 3 + 


8g 


bof 


i 6 A 5/2 (i) ♦ J ♦ sfcA 2 * + f x 5 l 5/2 ♦ f A 5 ^ 

e"* 2 + 36AL 5 / 2 + »-Oa £ L 5 0 + A 3 L T ^ 2 0 2 + ^ a 5 L T//2 ) - 

e“erfc 0 (jL 2 + 36AL 5 / 2 0 + 4OA 2 L 3 0 2 + 20A 2 L 3 + ^ A 3 L 7//2 0 3 + 
l6A 3 L T / 2 0 


Gi. = 


erfc ♦ - e^erfc 0^1 + l6A 2 L0 2 + 8* 2 L + ^ A 3 L 3/2 0 3 + ^ A^lV + 
32A 4 l 2 0 2 + 8A 4 L 2 + 32A 3 L 3 ^ 2 0^ + -=jp ^48A 2 L0 + 64a 3 l 3/ ^ 2 0 2 + 64A 5 L 3 ^ 2 + 
32A 4 L 2 0 3 + SOA^L^e” 1 ^ 


^ = e “^ 2 J-(48AL 3//2 + 144A 2 L 2 0 + 128A 3 L 5 / 2 + 128A 3 L 5 / 2 0 2 + 32A 4 L 3 0 3 + 
d| 3l/« \ 


8OA^L 3 0 j - e^erf c 0 j + |l6AL 0 + 48A L 0 2 + 24A 
64A 3 L 5 ^ 2 0 + ^ A 4 lV + 32A^L 5 0 2 + 8 aV + L j 


2 l 2 + 128 a 3 l 5/2 0 5 + 


3 
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= e“^ L 3 / 2 + 64AL 5 / 2 + 112X 2 L 3 0 + 64A 3 L T ^ 2 0 2 + 64A 3 L 7 / 2 + 

fcn 2 

A L^0 3 + y X L fij - e^erfc 0^9L 2 + 64AL 3 / 2 0 + 112A 2 L^0 2 + 
56 a 2 L 3 + 64A 3 L 7 / 2 0 3 + 96x 3 L 7 / 2 0 + ^ aN^V* + 32A 1 + L + 0 2 + 8X + L^ 

= -erfc t - e^erfc 0(ll + 4 a/l0 + 16X 2 L0 2 + 8A 2 L + a 5 L 5 / 2 0 3 + 

64A 3 L 3 / 2 0 + 32X 4 L 2 0^ + 96X 4 ‘L 2 0 2 + 24A^L 2 + — A 5 L 3 / 2 0 5 + 

15 

i|^ A 5 L 5 / 2 0 3 + 327?l?/ 2 $j + JL^4A\/l + l6A 2 L0 + A 5 L 3 / 2 0 2 + 

^ A 3 L 3//2 + 32A^L 2 0 3 + 8OA J+ L 2 0 + -hPjptyb + 1|2 a 5 l 5/202 + 

|i a 3 !?/ 2 ^-* 2 

^5. = _e a erfc 0 i|6OAL 3//2 0 + 24OA 2 L 2 0 2 ’+ 120A 2 L 2 + + 

48oa 3 l 3 / 2 0 + i6oa 1+ l 3 0 j+ + 48 oa^l 5 0 2 + i20a\ 3 + ^ a 3 l 7 / 2 05 + 
128A 3 L 7//£ 0 3 + 96A 5 L 7 / 2 0 + 3Lj + ^=^60AL 3//2 + 24OA 2 L 2 0 + 320A 5 !?/ 2 

32OA 5 L 5 // 2 0 2 + l6OA^L 5 0 3 + 4 OOa\ 3 0 + ±|S TpL?/ 2 ^ + ^ A 3 L 7 / 2 0 2 + 
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= -e Q erfc 0 (llL 2 + lOOAlP/ 2 0 + 24OA 2 L 3 0 2 + 120A 2 L 5 + A 3 L?/ 2 0 5 + 

So) \ 

32OA 3 L^/ 2 0 + %ik + 224A^L^0 2 + ^A^ + A 5 L 9 / 2 05 + 

A 5 L 9 / 2 0 3 + 32 A 5 L 9 / 2 0 ^ + -== ^JOOAL 5 / 2 + 720 A 2 L 5 0 + 64oa 5 iJ/ 2 0 2 + 

640A 5 L T / 2 + -224A if L 4 0 3 + 560 X ^0 + ^ A 5 L 9 / 2 0^ + A 5 L 9 / 2 0 2 + 

a 5 l 9/2 + 3L 3/2 ^-t 2 
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APPENDIX C 


NUMERICAL COMPUTATIONS EOR THE VARIOUS G n , 


dG, 




S, AND 


A 


SOn 

Si 


The c omp utations necessary to evaluate the equations of appendix B 
were carried out on a card-programed calculator. On this machine , the 
n umb er of significant figures in any computational step was limited to 
eight, a limi tation which resulted in a progressive diminishing of the 
significant figures in the answers as the value of 
given X and L. 


fact that each answer ^except for 


n increased at a 
The reason for the loss of significant figures is the 

^2.^ is the difference between two 


numb ers; thus, when these numbers become large and nearly equal, the 
resultant difference has few significant figures . 


It was also found that the complimentary error functions had to be 
computed to within a very small fraction of a percent of the correct 
value in order to obtain usable accuracy. The above-mentioned subtraction 
of two large but nearly equal numbers, one of which contains the compli- 
mentary error function .as a multiplier, is the reason for the particular 
accuracy required. The following expressions were used to compute the 
complimentary error function of 0 and \|r in the ranges of the argu- 
ments listed (0 is taken as the illustrative argument): 

erfc 0 = -^eL e - ^ 2 ^^ + a^T] 2 + a^T]^ + a^T}^ + a^r (0 0 < 2) 

(Cl) 


where 


1 

n 1 + O.32759H0 

a x = 0.22583685 

a 2 = -0.25212867 
a 5 = 1.2596951 
ail . = -1.2878225 
a^ = 0.9^064607 
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erfc 0 =0.000^0695202 - 0.00217828426 (l - 2.5e + 3-8333333e 2 - 3-9583333e 5 + 
2.80833336^ - 1. 2847222s 5 + 0. 24900794s 6 + 0.119667666^- 
0.1l490024e 8 + 0.036l75871e 9 ) (2 ^ 0 ^ 2.8) (C2) 


where 6 = 0 - 2.5 and 


erfc 0 = = 


-02 


1 l 1-3 1-3-5 A 1-3-5-7 1-3-5-7-9 


\fn 0 


- i + W- 

X yp 


x3 


x5 


1.3.5.7.9.11 1.3.5.7.9.11.15 


x 6 

76.331 104. 18) 

* f 1 


1 + - 
(2.8 < 0 < 00) 


(C3) 


p 

where X = 20 . 

Equation (Cl) is the approximation found on page 169 of reference 6. 
The errors Introduced by this approximation became unacceptable in this 
application above a value of the argument of 2. Equation (C3) was orig- 
inally thought to be useful in the range of arguments from 2 to <» since 

the product of the first bracketed expression and — 2— — represent the 

0 

first eight terms of the semiconvergent series for erfc 0 valid for 
large 0 and the constants of the second bracketed expression are chosen 
to give true answers at 0=2, 2.5, and 3- However, it was discovered 
that this equation was also unacceptable in .the range 2 < 0 < 2.8. A 
usable equation (eq. (C2)) for this range was found by employing the first 
10 terms of the Taylor expansion of erfc 0 about the point 0 = 2.5* 
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INSULATED PLATE AT VARIOUS VALUES OF L AND A 


NACA TN 4308 


1 
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9.6014 
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6.9730 
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5.2706 
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- 3-6790 

3.6776 
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-2.3698 

2.3690 

.76290 
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.85906 

- 1.4092 

1.4086 

.85901 
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1.1313 x 10- 10 
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6.9764 

5.2732 

3 . 6791* 

2.3702 

1.4093 

7.7360 x 10-5 
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INSULATED PLATE AT VARIOUS VALUES CF L AND X - Continued 


■ 

t - 0 

i = 1.0 


T 1 

& 

a& 

T 1 



T e 

T e 


T e 

T e 



Si 

&> 2 


io 2 

L - 0.1 

BBS 

0.03472 

-9.6528 X 10-2 

0.16877 

0.00078 

2.1252 X 10-2 

■ .2 

.06906 

-9.3094 

9.3568 X 10-2 

.02286 

9.1054 

.4 

.17142 

-8.2858 

8.0172 

.12964 

8.4214 

.6 

.31720 

-6.8280 
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r 13594 
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-7.3176 X 10-2 
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7.5704 x io* 2 
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4.5512 
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1.4 
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-2.9849 
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.91493 

-1.7014 

1.5938 
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1.7558 
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8.4292 x 10 -3 

•95044 

9.2860 X IO" 3 

2.0 

•97792 
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4.1364 

.97568 

4.5564 

2.5 

.99732 

-5.3660 x 10-^ 

5.0269 x 10-^ 

.99704 

5.5372 X 10-^ 

3.0 

.99980 

-4.0820 x 10-5 

3.8237 x 10-5 

.99978 

4.2120 X 10-5 

5-5 
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-1.9445 x 10 - 8 

I.8215 x IO -8 
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4.0 
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5.4292 x IO -8 
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5.9808 x io -8 
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0.1 

0.01326 

-9.6674 x 10*5 

9.8358 x IO* 5 

0.00834 

9.8848 x IO -3 

.2 

.04234 

-9.5766 

9.5460 

.03756 

9.5936 

A 

.15030 

-8.4970 

8.4698 

.14606 

8.5120 

.6 

.30388 

-6.9612 

6.93,90 

.30040 

6.9736 

.8 

.47864 

-5.2136 

5.1970 

.47604 

5.2228 

1.0 

.63218 

-3.6782 

3.6664 

.63034 

5.6846 

1.2 

.76278 

-2.3722 

2.3664 

.76160 

2.3764 

1.4 

.85873 

-1.4127 

1.4082 

.85802 

1.4152 

1.6 

.92232 

-7.7678 X 10*^ 

7.7430 x 10*^ 

.92193 

7.7816 X 10"^ 

1.8 

.96056 

-3.9402 

3.9314 

.96036 

3.9510 

2.0 

.98154 

-1.8459 

1.8400 

.98142 

1.8522 

2.5 

.99804 

-1.9628 x 10*5 

1.9566 x 10*5 

.99803 

1.9633 X 10-5 

5.0 

.99987 

-1.2659 x 10*6 

1.2618 x 10-6 

.99987 

1.2682 X IO* 6 

3.5, 

1.0 

-4.9586 X 10" 8 

4.9428 X IO* 8 

1.0 

4.9674 X IO -8 

4.0 

1.0 

-1.1805 x IO' 9 

1.1767 x io -9 

1.0 

1.1826 X 10*9 

L - 0.05 

0.1 

0.02484 

-4.8758 x io* 2 

6.1436 X 10 " 2 

0.00308 

3.5584 X IO" 2 

.2 

.05446 

-4.7277 

4.6533 

.03070 

4.7633 

.4 

■15972 

-4.2014 

4.1322 

.13860 

4.2361 

.6 

.30976 

-3-4512 

3.3944 

.29241 

3.4794 

.8 

.47593 

-2.6204 

2.5772 

.46276 

2.6420 


.65219 

-1.8390 

1.8088 

.62294 

1.8543 


.76140 

-1.1930 

1.1734 

.75540 

1.2029 


.85692 

-7.1540 x 10*5 

7.0362 X 10‘ 3 

.85333 

7.2128 X IO -3 

fill 

.92070 

-3.9648 

3.8996 

.91871 

3.9975 

1.8 

.95937 

-2.0313 

1.9979 

.95835 

2.0481 

2.0 

.98076 

-9.6190 x lXT^ 

9.4607 X 10*^ 

.98028 

9.6982 X 10-* 

2.5 

-99790 

-1.0522 

1.0349 

•99784 

1.0608 

3.0 

•99986 

-7-0375 x 10-6 

6.9217 x 10-6 

•99986 

7.0953 x IO- 6 

3-5 

.99999 

-2.8798 x 10*7 

2.8524 x 10-7 

.99999 

2.9036 x 10-7 

4.0 

1.0 

-7.2015 x 10 " 9 

7.0830 x IO* 9 

1.0 

7.2610 x io -9 
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TABLE II.- VALUES CF 


X 1 
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i = 0 AND 1.0 FOR PERFECTLY 


INSULATED PLATE AT VARIOUS VALUES CF L AND X - Continued. 
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S - 0 

I - 1.0 

T 1 

Te 

a 

_£e 

Je 

StD 2 

E 1 

T e 

aS, 

T e 

io 2 

t-« 

N 

O 

Ch 

0.1 

0.05890 

-0.28233 

0.84238 

0 

0.00103- 

.2 

.11248 

-.26626 

•38414 

.00586 

.13318 

.4 

.21628 

-23512 

.21551 

.09634 

.24388 

.6 

.34654 

-.19604 

■17791 

.24622 

.20522 

.8 

.49315 

-.15205 

•13799 

.41534 

.15918 

1.0 

.63441 

-.10968 

9.9534 X 10- 2 

.57829 

.11482 

1.2 

.75476 

-7.3572 x 10* 2 

6.6768 

.71712 

7.7016 X 10 “2 

1.4 

.84703 

-4.5891 

4.1648 

.82355 

4.8041 

1.6 

.91126 

-2.6622 

2.4160 

.89763 

2.T870 

1.8 

.95212 

-1.4363 

1.3035 

.94477 

1.5036 

2.0 

.97598 

-7.2066 X 10-3 

6.5400 x 10-3 

.97229 

7-5444 x 10-3 

2.5 

.99888 

-9.3516 x 10-4 

8.4870 x 10-4 

.99640 

9.7896 x 10-4 

5-0 

-99974 

-7.7094 x 10* 5 

6.9966 x 10 -5 

.99970 

8.0706 X 10“ 5 

3-5 

■99999 

-4.0385 x 10- 6 

3.6650 x 10-6 

.99998 

4.2277 x 10-6 

4.0 

1.0 

-1-3423 x 10-7 

1.2182 x 10-7 

1.0 

1.4052 x 10-7 

L - 0.4 

0.1 

0.06752 

-0.37299 

1.2780 

0 

0.00017 

.2 

.12808 

-34877 

.57422 

.00288 

.10955 

.4 

.23690 

-.30524 

.27893 

.08192 

: 31266 

.6 

.36087 

-.25565 

.22498 

.22917 

.27124 

.8 

-50043 

-.19982 

.17582 

•39748 

.21205 

B3 

.63605 

-.14558 

.12809 

.56105 

'.15448 

m 

.75287 

-9.8848 x 10- 2 

a. 6968 x io-2 

.70194 

.10490 

m 

.84362 

-6.2554 

5.5036 

.81138 

6.638O X lO- 2 


.90775' 

-3-6899 

3.2465 

.88874 

3.9154 

1.8 

.94929 

-2.0286 

1.7848 

•93883 

2.1527 

2.0 

• 97402 

-1.0394 

9.1448 x 10" 3 

.96866 

1.1030 

2.5 

•99641 

-1.4357 x 10-3 

1.2631 

.99567 

1.5235 X 10-3 

3.0 

•99968 

-1.2773 x 10-4 

1.1238 x 10-4 

•99961 

1.3554 x 10-4 

3.5 

•99998 

-7.3174 x 10-6 

6.4381 x 10-6 

•99998 

7.7650 x 10-6 

4.0 

1.0 

-2.7017 x 10-7 

2.3770 x 10-7 

1.0 

2.8670 x 10-7 






70 


NACA TN 4308 


c£i 

TABLE II.- VALUES OF ^1, -k, AND -k AT £ 


6 | 




INSULATED PLATE AT VARIOUS VALUES OF L 


= 0 AND 1.0 FOR PERFECTLY 
AND A - Concluded 


A 

£ = 0 

| = 1.0 

T 1 

T e 

& 

8| 

& 

Jk 

8ii 2 

T 1 

T e 

T, 

d-i 

L = 0.5 

0.1 

0.07504 

-0.46248 

1.7633 

0.00003 

0 

.2 

. 14151 

-.42924 

.78278 

.00152 

.08179 

.4 

.25600 

-.37200 

.34551 

.06922 

.36890 

.6 

.37499 

-.31250 

.26710 

.21303 

.33547 

.8 

.50787 

-.24606 

.21002 

.38028 

.26447 

1.0 

.63796 

-.18102 

.15450 

.54409 

.19456 


.75139 

-.12430 

.10610 

.68693 

. 13361 


.84049 

-7.9755 x 10 " 2 

6.8072 X 10 -2 

•79914 

8.5719 X 10“ 2 


.90442 

-4.7792 

4.0791 

.87964 

5.1365 

1.8 

.94650 

-2.6748 

2.2830 

.93264 

2.8748 

2.0 

.97204 

-1.3982 

1.1933 

.96479 

1.5027 

2-5 

.99510 

-2.0490 x 10 " 5 

1.7488 x 10 

.99484 

2.2022 X 10 "3 

3-0 

.99961 

-1.9596 X 10-^ 

1.6725 x 10 

.99951 

2.1061 X 10 

3.5 

•99998 

-1.2226 X 10“ 5 

1.0436 X 10“ 5 

•99997 

1.3141 X 10 "5 

4.0 

1.0 

-4.9836 X 10 -T' 

4.2536 x 10 " 7 

1.0 

5.3563 x io" 7 
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Figure 2.- Plot of the nondimens ional parameter L as a function of 
skin material and thickness for h = 0.1 Btu/sq ft-°F-sec. 
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Figure J.- Variation of the parameter p with the parameter (pck), for various outer skins 
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hacking-material curve with the abscissa for the outer skin. If the outer skin and backing 


material are interchanged, the sign of p is reversed. 
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Figure 4.- Continued. 
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Figure 4.- Concluded. 




Temperature ratio 


1.2 



(a) L = 0.001. 

Figure 5*- Ratio of T-Jl^ a t 5=0 for a finite number of terms in equation (25) to 
from equation (15) plotted against the parameter A. 0=1. 
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(c) L = 0.5. 
Figure 5.- Concluded. 
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Figure 6.- Temperature distribution at the outer face plotted against time for composite slabs 
having stainless -steel skins suddenly exposed to aerodynamic heating at t = 0. The skin 
thicknesses Eire 0.060 and 0.030 inch and the insulating backing materials are mahogany, balsa, 
and a perfect insulator, h = 0.056 Btu/sq. ft-°F-sec. 
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Figure 7*- Ratio of indicated heat-transfer coefficient to true heat-transfer coefficient plotted 
against time for a composite slab having a 0.050-inch-thick stainless-steel skin and various 
insulating hacking materials, h = 0.056 Btu/sq. ft-°F-sec. 
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Figure 8.- Ratio of indicated heat-transfer coefficient at 5=0 to true heat-transfer coeffi- 
cient plotted against time for composite slabs having stainless- steel skins of thickness l 
and various insulating backing materials, h = 0.056 Btu/sq ft-°F-eec. 
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Figure 9 .- Ratio of indicated heat- transfer coefficient to true heat- 
transfer coefficient plotted against ratio of local temperature to 
equilibrium temperature for composite slab having a 0 . 030 -inch-thick 
stainless-steel outer skin and various insulating backing materials, 
h = 0.056 Btu/sq. ft-°F-sec. 
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